Effects of Isotope Substitution on Local Heating and Inelastic Current in Hydrogen 

Molecular Junctions 
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Using first principle approaches, we investigate the inelastic features in the hydrogen molecular 
junction. We observe that local heating and inelastic current have significant isotope-substitution 
effects. The junction instability is also relevant to the isotope substitution. We predict that the HD 
junction has the smallest breakdown voltage compared with the H2 and D2 junction in the optimized 
geometry. The selection rule for modes that significantly contribute to the inelastic effects is related 
to the component of vibration along the direction of electron transport. 



Building electronic circuits from molecules is an in- 
spiring idea. [1] In the past decade, motivated by this 
revolutionary idea, intensive theoretical and experimen- 
tal studies have been conducted to explore the charge 
transport in a single molecule junction. Single molecule 
junction is the basic building block for electronic compo- 
nents in molecular electronics. Understanding the trans- 
port properties of the basic building block is of key im- 
portance in the design of a new form of electronic circuit 
at the atomic and molecular levels. [3, 0, 0, [1]. 

The hydrogen junction is an ideal testbed to compare 
theories and experiments. Despite the simplicity of the 
hydrogen junction, the data reported for the conductance 
varying from 0.2 to 1.0 Go have provoked controversy. 
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FIG. 1: (color online) Schematic diagrams for the systems 
investigated: (a) the optimized geometry of the H2 junction. 
The Pt-H-H-Pt atoms are sandwiched between two bulk Pt 
electrodes kept at a certain external bias with the chemical po- 
tentials of the left/right electrodes are shifted, i.e., the source 
drain bias is V sd=(Eifr-Efl) /e. The Pt electrodes are de- 
scribed as Jellium model with the interior electron density of 
Pt electrodes is taken close to the value of metallic Pt (r s ~ 3). 
The H-H, Pt-H, and Pt- Jellium distance are 1.615, 4.051, and 
2.0 a.u., respectively; (b) a specific geometry of the H2 junc- 
tion. The hydrogen molecule rotates an angle of 15° due to 
asymmetric contacts. 



H 0, 0, 0, E3, EH The problem of inelastic electron- 
vibration interactions remains unsettled. For example, 
the features of current-voltage characteristics at a bias 
around 65 mV are attributed to the electron-vibration 
interactions. However, a theoretical group has re- 

ported calculations to support the idea that the static 
atomic structure can still cause a pronounced resonance 
at this bias. [6] The mode identification and the selec- 
tion rule for modes that significantly contributed to the 
inelastic electron tunneling spectroscopy (IETS) have not 
reached unanimity. The identification of the spike in 
dl/dV at the Pt-H2 system as the transverse vibration 
mode is still controversial. [l3[ 

To clarify the fundamental inelastic electron transport 
properties, we performed comparative studies on inelas- 
tic profiles among the H2, HD, and D2 junction. The 
simplicity of the hydrogen junction allows the inelastic 
profiles to provide clear information for understanding 
how the vibration modes contribute to inelastic effects in 
the nano scale junction. Owing to the small mass of the 
hydrogen atom, isotope substitution has a great influ- 
ence on inelastic effects. While the isotope substitution 
does not affect the contact geometry and the bonding 
strength, it can serve as a controllable factor that is sen- 
sitive enough to explore the fundamental properties of 
inelastic effects in the nano scale junction. Moreover, 
the vibrations excited before the breakdown voltage of 
the hydrogen junction are acoustic-like modes connected 
to contact geometry and bonding strength. This is in 
sharp contrast to the IETS and HREELS in the alkyle 
junction where all observed modes in experiments are 
optical-like modes more relevant to internal vibrations of 
the molecule. [H, [3, E3, H H, H EH As such, inelastic 
profiles in the hydrogen junction can provide more infor- 
mation on contact geometry and bonding strength in the 
molecular junctions. Although heating is an important 
effect relevant to the stability and performance of device, 
very few attempts have been made that addressed the 
nano junction. [14[ 

In our results, we observed interesting inelastic fea- 
tures in the hydrogen junction. In the optimized ge- 
ometry (Fig. Hl[a)), The electron- vibration interaction is 
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FIG. 2: (color online) Schematic diagrams for the vibrational 
modes of the H2 molecule. Six vibrational modes are shown 
in (a) for the modes 1 and 2 associated with a rigid motion 
of the H2 molecule perpendicular to the line connecting two 
Pt atoms; (b) for the mode 3 associated with a rigid motion 
of the H2 molecule along the line connecting Pt atoms; (c) 
for the modes 4 and 5 associated with a transverse motions 
of the H atoms; (d) for the mode 6 associated with an inter- 
nal vibration of the H2 along the line connecting Pt atoms. 
The energies and vibration characteristics of these modes are 
shown in Table I for the linear junction with a geometry de- 
picted in Fig. [TJa). 



TABLE I: Vibration Modes depicted in Fig. [2] 

modes lvh 2 uhd (£^) vd 2 (^0 features 

1 36.8 30.0 (1.22) 26.1 (1.41) transverse 

2 36.8 30.0 (1.22) 26.1 (1.41) transverse 

3 64.8 52.9 (1.22) 45.9 (1.41) longitudinal 

4 115.9 101.1 (1.15) 82.1 (1.41) transverse 

5 115.9 101.1 (1.15) 82.1 (1.41) transverse 

6 314.9 273.4 (1.15) 222.7 (1.41) longitudinal 
Unit of energies in meV 



elastic current and heating effects in the linear hydro- 
gen junction (Fig. [2(a)). We study these effects by the 
first order perturbation theory based on the second quan- 
tized formalism where the wave functions are obtained 
from density functional theory (DFT). [22[ The scatter- 
ing wave functions of the whole system are calculated 
by solving the Lippmann-Schwinger equation iteratively 
until the self-consistency is obtained. [24] The exchange 
and correlation energy are calculated at the level of the 
local density approximation. [25] These wave functions 
are applied to calculate current and heating. The heating 
phenomenon in nano junctions has been measured experi- 
mentally in the alkane junctions in C-AFM, and the data 
are in reasonable agreement with our theoretical calcula- 
tions, mm 

At first glance, the conductance of the hydrogen 
molecule is expected to be small because of a closed shell 
configuration. However, the hybridization of states be- 
tween hydrogen and Pt atoms can create new channel for 
current. [7[ As a result, the magnitude of conductance, 
irrelevant to the isotope substitution, can be quite large 
around 1.0 Goin the hydrogen junction. This feature is 
verified by the analysis of partial DOS. In our results, the 
current- volt age characteristic in the Pt-H-H-Pt junction 
is quite linear, with a conductance around G « 0.9 Go- 
We find that the S-wave state, mostly contributed by the 
hydrogen atoms, has negligible contribution to the con- 
ducting state. This feature supports the argument that 
the origin of charge transport comes mainly from the H- 
Pt hybridized states. Q 

We briefly introduce the theory of the electron- 
vibration interactions. By applying the field operator 
with the wave functions obtained from the static DFT, 
the electron- vibration interaction has a second quantized 
form, [H 



suppressed by the mass of the molecule. The electron- 
vibration interaction is found to be enhanced by asym- 
metric molecule, which leads to the most prominent heat- 
ing and inelastic correction to the current in the HD 
junction. Since the current-induced force, irrelevant to 
the isotope substitution, is found to be small, we pre- 
dict that the HD junction is the most unstable owing to 
local heating as the major breakdown mechanism. This 
feature can be verified by experiments. In a possible ro- 
tated geometry as shown in Fig. [2(b), we observe that all 
modes significantly contribute to inelastic profiles. It fur- 
ther confirms that the selection rule of important modes 
in the inelastic profiles is connected to the component of 
molecule vibration along the direction of charge trans- 
port. Thus, inelastic profile can a sensitive tool to probe 
the contact geometry and bonding strength in the nano 
scale junction. 

We will begin our discussion by considering the in- 
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where Mi is the mass of the i — th atom, Ai^j^ is a canon- 
ical transformation between normal coordinates and 
Cartesian coordinates satisfying V\ A^ j v A- •/ / 



^•,j /(S ^y 5 ^ % Ei°E2 * s coupling constant related to the 
electron and vibration of atoms, which can be calculated 
from the pseudopotential of atoms and wavefunctions; 
a^' R ^ is the annihilation operator for an electron with 
energy E incident from the left or right electrode, respec- 
tively; fajjj V represents the energy of the jV-th normal 
mode and bj V is the annihilation operator corresponding 
to that mode. The strength of electron-vibration cou- 
pling is affected by the isotope substitution related to 
the mass of the molecule shown in Eq. (pQ). 
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FIG. 3: (color online) Differential conductance (left axis) and 
absolute value of d 2 I/dV 2 (right axis) as a function of bias 
for the H2 (black solid line), HD (red dotted-dashed line), 
and D2 (blue dashed line) molecule sandwiched between the 
optimized linear junction (upper panel) and the junction of 
with asymmetric contacts (lower panel). The schematics show 
the noticeable modes contributed to the inelastic current in 
the H2 junction. 



As the first step in our analysis, we focus on inelas- 
tic profiles in an optimized linear geometry as shown in 
Fig. UK a). Two heavy Pt electrodes, which are inert ow- 
ing to heavy masses, are simulated as two Pt atoms with 
infinite mass for simplicity. The geometry is relaxed in 
gauusian03 with the total energy of the system around 
—380.15 Hartree. The geometry is stable and requires 
0.53eV to take the H2 molecule away from the junction. 
In Table H we list the energies of all the modes, while the 
schematics of vibrations for the H2 molecule are shown 
in Fig. [21 Mode 3 is the one responsible for significant 
inelastic current and heating because of its longitudinal 
nature. Mode 6 is associated with a longitudinal internal 
vibration. The energy of this mode is too large to be ex- 
cited before the breakdown bias around 200 mV found in 
experiments, d, Q The rest of the modes have insignif- 
icant inelastic effects owing to their transverse natures. 

In the linear hydrogen junction, the position of the 
discontinuities in the differential conductance is affected 
by the isotope substitution as shown in the upper panel 
of Fig. [3l In the H2 junction, an abrupt change in 
the differential conductance appears at a bias around 
65 mV in the d 2 I/dV 2 spectrum. This peak is associ- 
ated with mode 3 which has longitudinal character. This 
mode is associated with motion of the molecule as a rigid 
body. Accordingly, the energies are perfectly scaled by 
^h 2 / u °hd v thh d I mH 2 ~ 1.22. The rest modes with 
transverse characters are too weak to be resolved. 

The magnitude of electron- vibration interaction is also 
affected by the isotope substitution. Let us compare the 
magnitude of discontinuities in the differential conduc- 
tance in the linear H2 and D2 junctions. These two junc- 
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FIG. 4: (color online) Local temperature as a function of bias 
for the H2 (black solid line), HD (red dotted-dashed line), 
and D 2 (blue dashed line) molecule sandwiched between the 
optimized linear junction (upper panel) and the junction with 
asymmetric contacts (lower panel). The schematics show the 
noticeable modes contributed to heating in the H2 junction. 



tions have the same mechanical structures except for dif- 
ferent masses. The electron- vibration coupling is related 
to the mass of ions by 1/y/Mi. Thus, the H2 junction has 
a larger discontinuity in the differential conductance be- 
cause of smaller mass. In the HD junction, the electron- 
vibration interaction is enhanced by asymmetric mass 
distribution. As a result, the HD junction has the largest 
change in the differential conductance. 

Local temperatures in the linear H2, HD, and D2 junc- 
tions are further investigated to elucidate the role of iso- 
tope substitution. Two major processes lead to an equi- 
librium local temperature in a nano junction. One is due 
to the electron-vibration interaction that occurs in the 
atomic region of the junction; the other is due to the dis- 
sipation of heat energy to the bulk electrodes via thermal 
transport. [27[ We assume that the energy generated in 
the atomic region via inelastic electron-vibration scat- 
tering and the energy dissipated to the electrodes via 
thermal current finally reach equilibrium, such that a 
well-defined local temperature can be calculated in the 
atomic region. We evaluate the probability of the energy 
exchange via electron-vibration interactions by the Fermi 
golden rule for the four first-order scattering mechanisms. 
These mechanisms correspond to electrons incident from 
the right or left electrode that relax (cool) or excite (heat) 
the energy level of the normal mode vibrations in the 
atomic region. The total thermal power generated in the 
junction can be calculated as the sum of all vibrational 
modes of the above four scattering mechanisms. Detailed 
theories can be found in Ref. [22[. 

Heating induced by the electron-electron interaction is 
neglected in present work since the contribution might 
be small compared with the electron-vibration interac- 
tion at small biases. 3, IK In the upper panel of Fig. [H 
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we plot the equilibrium temperature of the molecule as a 
function of applied voltage in the linear H2, HD, and D2 
junctions. Weak heating phenomena occur at small bi- 
ases corresponding to modes 1 and 2. These modes have 
very small contribution to heating due to their transverse 
vibration characters. Major heating profiles are observed 
at onset biases around 64.8 mV, 52.8 mV, and 45.9 mV 
for the H2, HD and D2 junction, respectively. Those 
heating profiles are related to mode 3 with longitudinal 
vibration character. 

The interplay between the onset bias and the strength 
of electron-vibration interaction caused by the isotope 
substitution leads to interesting heating profiles. For 
symmetric molecules in the linear junction (Fig. QJa)), 
the magnitude of electron- vibration coupling is enhanced 
by the smaller mass (the H2 junction) which leads to a 
larger heating rate, whereas the molecule with a larger 
mass (the D2 junction) has a smaller onset bias. Con- 
sequently, a crossing in the equilibrium temperatures for 
the H2 and D2 junctions occur at 80 mV as shown in the 
upper panel of Fig.JH For the biases smaller than 80 mV, 
the local temperature in the D2 junction is larger due to 
a smaller onset bias. For the biases larger than 80 mV, 
the local temperature in the H2 junction is larger due 
to the larger electron-vibration coupling. For the HD 
junction, the asymmetric mass distribution leads to the 
enhancement of electron-vibration interaction. It accord- 
ingly induces the largest heating rate after the longitudi- 
nal normal mode is excited. 

The junction stability is also relevant to isotope sub- 
stitution. In atomic-sized Al wires, the current-induced 
force is the major mechanism that causes mechanical 
breaking, whereas in alkane junctions heating is the dom- 
inant breaking mechanism. [22I, 0] In our results, the 
current-induced forces on the hydrogen atoms have a 
small value of 0.38 nF at 200 mV. We thus expect that 
the major breakdown mechanism will be caused by local 
heating rather than the current-induced force. Because 
local heating is susceptible to effects of isotope substitu- 
tion in contrast to the current-induced force where the 
isotope substitution plays no role, in the optimized ge- 
ometry (Fig. QJa)), we predict that compared with the 
H2 and D2 junctions, the HD junction has the smallest 
breakdown voltage caused by the most prominent heat- 
ing. 

As the detailed configurations in the contact region 
are unknown in the real experiment, the optimization 
of contact is not our present concern. We will concen- 
trate on how the orientation of the molecule affects the 
inelastic profiles. We may now proceed to consider a 
specific configuration of the hydrogen junction where the 
hydrogen molecule rotated at an angle of 15° as shown in 
Fig. DJb). The modes with transverse characters in the 
linear junction turn out to have large component of vibra- 
tion along the direction of current in rotated geometry. 
The responses of inelastic current and local heating to 



the biases are shown in the lower panels of Figs. [3] and [4j 
respectively. In contrast to the optimized case (Fig.[]Ja)) 
where only mode 3 had significant contribution to inelas- 
tic effects, all modes that can be excited cause prominent 
features in the inelastic profiles. These features further 
illustrate the selection rule for modes that significantly 
contribute to the inelastic effects related to the compo- 
nent of vibration along the direction of electron trans- 
port. 

In conclusion, we have calculated various transport 
properties in H2, HD, and D2 molecular junctions in lin- 
ear and rotated geometries. In the linear geometry, the 
I-V characteristics is linear with conductance G ~ 0.9 Go- 
We investigate how the isotope substitution can influence 
the inelastic current and local heating. The interplay be- 
tween vibration energies and the electron- vibration cou- 
pling strength leads to interesting heating profile. We 
predict that the HD junction has the smallest breakdown 
voltage due to the enhancement of electron- vibration in- 
teraction. This feature is measurable and could give an 
insight into the mechanisms responsible for the system 
breakdown. We also investigate the inelastic effects in ro- 
tated geometry. The inelastic features further illustrate 
the selection rule for modes that significantly contribute 
to the inelastic effects related to the component of vibra- 
tion along the direction of electron transport. 
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